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ABSTRACT 


This  report  presents  a  user 1 s  manual  for  a  computer 
program,  UCIN-CABLE:  Underwater  Cables.  The  program 
is  designed  to  study  the  dynamics  of  underwater  cables. 
The  cable  is  modeled  as  a  series  of  riaid  cylinders 
connected  end-to-end  by  ball  and  socket  ioints.  There 
may  be  one  or  many  branches,  and  no  closed  loops  are 
formed.  The  effects  of  normal  fluid  draq,  tangential 
fluid  drag,  added  mass,  and  buoyancy  are  included. 

This  manual  provides  detailed  instructions  for 
building  a  set  of  input  data  and  interpreting  the 
output  of  the  program.  _ 
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I.  INTRODUCTION 


This  user's  manual  is  intended  to  accompany  the  computer 
program  UCIN-CABLE:  Underwater  Cables.  Subsequently,  the 
computer  program  is  intended  to  provide  numerical  data  on  the 
dynamics  of  underwater  cables.  The  program  models  a  flexible 
cable  by  a  series  of  rigid  cylinders  connected  end-to-end  by 
ball  and  socket  joints  [1,2]*.  The  size,  shape,  and  mass  of  each 
of  the  segments  is  arbitrary.  (See  Figure  1.)  The  program  also 
allows  for  spherical  anchors  and/or  towed  spheres.  (These  bodies 
may  only  be  at  the  end  of  the  branches.)  The  fluid  effects  on 
the  cable  (normal  fluid  drag,  tangential  fluid  drag,  added  mass, 
and  buoyancy  [2])  and  on  the  appended  spheres  (normal  fluid  dracr, 
added  mass,  and  buoyancy)  may  all  be  included  individually  as 
desired.  Finally,  the  program  is  expected  to  be  effective  in 
studying  the  nonlinear,  three-dimensional,  dynamical  behavior  of 
long,  heavy,  towing  and  hoisting  cables. 

The  dynamic  analysis  of  the  cable  model  is  a  specialization 
of  a  recently  developed  analysis  for  general  "chain"  systems 
[3-6].  (A  general  chain  system  is  a  set  of  rigid  bodies  arbitrarily 
assembled  and  connected  provided  that  adjoining  bodies  have  at 
least  one  point  in  common  and  that  no  closed  Ioods  are  formed.) 

In  fact,  except  for  the  fluid  force  modeling,  the  fundamental 
algorithms  of  UCIN-CABLE  are  the  same  as  for  the  previously  developed 
program  UCIN-SUPER  [7] .  The  dynamical  formulations  and  fundamental 

*  Numbers  in  brackets  refer  to  References  at  the  end  of  the  Report. 


algorithms  in  that  program  are  based  on  Lagrange's  form  of 
d'Alembert's  principle  or  Kane's  equations  [8-11]  and  other  vector- 
tensor-matrix  methods.  The  fluid  effects  were  simply  "added"  to 
this  formulation  to  form  UCIN-CABLE.  A  detailed  description  of 
the  dynamical  principles  and  methods  and  the  fluid  force  model 
appears  in  references  [1-6]. 

The  attractive  features  of  UCIN-CABLE 's  cable  model  are  as 
follows:  1)  the  arbitrary  dimensions  and  physical  parameters  of 
the  segments;  2)  the  arbitrary  specification  of  externally  applied 
forces;  3)  the  use  of  relative  orientation  annles  between  the  links 
to  define  the  system  configuration;  and  4)  the  use  of  Lagrange's 
form  of  d'Alembert's  principle  to  develop  the  governing  ecjuations 
of  motion.  Lagrange's  form  of  d'Alembert's  principle  has  been 
shown,  for  large  systems,  to  possess  the  advantages  of  both  La grange's 
equations  and  Newton's  law,  but  without  the  corresponding  disadvan¬ 
tages.  That  is,  the  principle  provides  for  automatic  elimination 
of  unknown  "non-working"  internal  constraint  forces  without  intro¬ 
ducing  tedious  differentiation  of  scalar  energy  functions.  Finally, 
the  governing  dynamical  equations  are  integrated  numerically  using 
a  fourth  order  Runge-Kutta  method. 

The  user's  manual  is  divided  into  four  sections  with  the  first 
describing  the  scope  and  range  of  applicability  of  the  program. 

The  second  section  describes  the  physical  data  and  governing 
variables  of  the  finite  segment  model,  and  the  data  for  the  surround¬ 
ing  water,  needed  to  run  the  program.  The  third  section  provides 
detailed  (card  image)  instructions  for  the  input  data  and  the  fourth 
section  describes  the  output  data. 
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II.  SCOPE  AND  RANGE  OF  APPLICATION 


The  program  is  designed  to  numerically  perform  the  following 
dynamical  analyses: 

1.  Given  the  physical  data  (sizes,  masses,  inertias,  and 
assemblage)  of  each  body  in  the  system,  oiven  the  forces 
applied  to  reference  point  of  body  1,  given  the  external 
forces  and  torques  applied  to  each  of  the  bodies,  and 
given  information  about  the  surrounding  fluid  (ambient 
velocity,  density,  viscosity,  surface  height,  etc.) 

the  kinematics  (position,  velocity,  acceleration,  angular 
velocity,  and  angular  acceleration)  of  each  body  is 
determined.  (Since  the  cable  is  assumed  perfectly 
flexible,  the  connecting  ball  and  socket  joints  are  free, 
and  hence  the  torques  between  adjoining  bodies  are  zero.) 

2.  Given  the  physical  data,  given  the  kinematics  of  the 
reference  point  of  body  1,  given  the  external  forces 
and  torques  applied  to  each  of  the  bodies,  and  given 
information  about  the  surrounding  fluid,  the  required 
driving  forces  for  the  reference  point  of  body  1  and  the 
kinematics  of  each  body  is  determined. 
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PHYSICAL  AND  FLUID  DATA  AND  GOVERNING  VARIABLES 


This  section  describes  the  physical  data  characterizing  the 
finite  segment  model  and  the  data  for  the  surrounding  fluid, 
needed  as  input  for  the  program.  The  variables  (coordinates) 
defining  the  kinematics  are  also  discussed  and  defined. 


1. 


Chain  Assemblage  -  Body  Connection  Array. 

Consider  again  the  chain  system  shown  in  Figure  1.  Let 
the  8  bodies  of  the  chain  be  labeled  B^  (i*l,...,8) 
with  B^ ,  the  reference  body,  being  chosen  as  the  segment 
at  the  beginning  of  the  cable.  The  other  bodies  are 
labeled  outward  along  the  "branches"  and  clockwise  as 
shown  in  Figure  2.  This  assemblage  is  entered  into  the 
program  through  an  array  called  JFV^  which  lists  for  each 
Bi  (i=l,...,N)  (for  N  bodies)  the  adjacent  lower  numbered 
body.  For  example,  for  the  chain  system  in  Figure  2.,  the 
JFV  array  is: 


JFV  =  (0,  1,  2,  3,  4,  2,  6,  7) 


where  0  (zero)  refers  to  an  inertial  reference  frame  R. 


Masses . 

The  mass  of  each  body  must  be  known  for  input  data.  The 
usual  units  are  slugs  or  kilograms  although  any  system 
of  units  may  be  used. 


.»***. -v*.*^fc  -’Safilfc  -*«  *•»  »*%*•  v  *4i' 


******  M**&**z4i^W«*m  A*M*  : 


rtutti 


3.  Inertia  Matrices. 

The  inertia  matrix  (moments  and  products  of  inertia) 
for  each  body  referred  to  mutually  perpendicular  body 
fixed,  coordinate  axes  passing  through  the  mass  center 

G^,  must  be  known  for  each  body  (i=l,...,N).  The 

2  2  2 
usual  units  are  slug-m  ,  slug-ft  ,  or  kg-m  . 

4 .  Mass  Center  Locations  --  R- Vectors. 

Let  0^  be  an  arbitrarily  chosen  reference  point  of 
Let  CK  be  a  reference  point  of  B^  (i=l,...,N)  and  let 
0^  be  the  common  point  of  B^  and  its  adjacent  lower 
numbered  body.  (See  Figure  3.,  for  example.)  Let 
be  the  mass  center  of  B^ .  Then  the  vector  ri  from  CK 
to  G^  locates  the  mass  center,  r^  is  fixed  in  Bi  and 
its  components  (required  as  input)  are  referred  to  the 
body-fixed  coordinate  axes  of  B^. 

5.  Reference  Point  Locations  —  £  -Vectors . 

Let  be  the  vector  from  the  reference  point  of  the 
adjacent  lower  numbered  body  of  B^  to  the  reference  point 
of  B^ .  Then  is  fixed  in  the  adjacent  lower  numbered 
body  and  its  components  (required  as  input)  are  referred 
to  the  body-fixed  coordinate  axes  of  the  adjacent  lower 
numbered  body  of  B^ .  (See  Figure  3.) 

6.  Variable  Definition. 

If  the  chain  system  contains  N  bodies,  3N+3  variables 
(coordinates)  are  required  to  define  the  position  and 
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configuration  of  the  system.  This  can  be  seen  by 
observing  that  three  variables  are  needed  to  locate 
0^  relative  to  the  inertial  frame  R;  three  are  needed 
to  define  the  orientation  of  in  R;  and  three  are  needed 
to  define  the  orientation  of  each  of  the  remaining  N-l 
bodies  relative  to  their  adjacent,  lower  numbered  bodies. 

Let  these  3N+3  variables  be  called  ( j=l, . . . , 3N+3) 
and  consider  them  in  succession  in  sets  of  three  each. 

Then  each  set  except  the  first  may  be  considered  as  three 
orientation  angles  cu  ,  8^  and  y^  of  relative  to  its 
adjacent,  lower  numbered  body.  These  orientation  angles 
are  discussed  in  detail  in  the  Appendix,  but  they  are 
basically  angles  obtained  by  three  successive  dextral 
rotations  of  B.  relative  to  its  adjacent,  lower  numbered 
body  about  its  body-fixed  coordinate  axes. 

7.  Specified  Motion  of  the  Bodies  of  the  System. 

If  the  motion  of  some  of  the  bodies  of  the  chain  system 
is  known,  the  corresponding  variables  (coordinates)  may 
be  specified  through  functions  of  time,  etc.  in  SUBROUTINE 
START  of  the  program. 

In  many  cases  (e.g.,  finite  segment  model  of  a  towed  cable), 
it  is  desired  to  specify  the  acceleration  of  some  of  the  bodies. 
Sometimes  the  desired  acceleration  is  known  simply  in  terms  of  a 
curve  or  data  points  rather  than  through  elementary  functions.  To 
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accomodate  this  situation,  the  program  will  "read  in"  the  second 
derivative  of  the  variables  at  selected  time  intervals.  The 
program  has  the  capability  of  including  25  points  of  an  acceleration 


curve . 


To  be  more  specific,  consider  the  acceleration  curve  and  the 
straight  line  approximation  shown  in  Figure  4.  The  acceleration, 
velocity,  and  displacement  during  the  ith  time  interval  are  then: 


a  ■  ‘  1  (t  *  V 


a.,  -a.  (t-t.) 

v  -  v •  +  a  (t  -  t.)  +  (  1  )  - 

ii  l  ci+1  ti 


d  = 


d.  +  v.(t  -  ti)  +  &i 


(t  -  ti)  .  ,  ai+l  -  ai 
ti+l  ~  ti 


+  ( 


where  a^  v^,  and  ti  are  the  acceleration,  velocity,  displacement, 
and  time  at  the  beginning  of  the  ith  interval.  The  motion  of  0^ 
the  reference  point  of  is  thus  determined  when  the  a.^  are 
specified  and  when,  v^,  and  d^,  the  initial  velocity  and  displacement 
are  given. 


Acceleration 


IV.  INPUT  DATA 


The  following  paragraphs  describe  the  card  image  requirements 
for  the  input  data.  An  example  listing  is  provided  in  the  next 
section. 

1.  START  OF  DATA  Card.  The  first  card  must  read  "*START  OF 
DATA*".  (See  example  input.) 

2.  Number  of  Bodies.  The  next  card  contains  the  number  of 
bodies  in  the  chain  system  NB,  read  in  15  format. 

3.  Mass  Center  Labels.  The  program  provides  an  option  of 
labeling  or  naming  the  mass  centers  of  the  bodies  of  the 
system.  Hence,  the  next  card  contains  the  number  of  mass 
center  labels  (names)  to  be  assigned,  read  in  15  format. 

If  no  labels  are  to  be  read,  a  blank  card  is  used. 

If  the  number  read  is  >.  1,  the  next  series  of  cards  then 
contain  the  mass  center  labels  (or  names)  (5  per  card) 
read  in  5A16  format  until  the  total  number  read  in  3 
is  exhausted. 

4.  Joint  Labels.  The  program  also  provides  an  option  of 
labeling  or  naming  the  connecting  joints  between  the 
bodies.  Hence,  the  next  card  contains  the  number  of 
joint  labels  (or  names)  to  be  assigned,  read  in  15  format. 
If  no  labels  are  to  be  read,  a  blank  card  is  used. 
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If  the  number  read  >_  1,  the  next  series  of  cards  then 
contain  the  joint  labels  (or  names)  (5  per  card)  read 
in  5A16  format  until  the  total  number  read  in  5.  is 
exhausted. 

5.  Output  Printing  Reference  Body.  Let  the  bodies  of  the 
chain  system  be  numbered  as  discussed  in  the  previous 
section.  Then  the  next  card  contains  the  body  number 
or  that  body  which  is  to  be  an  "output  printing" 
reference  body  of  the  system,  read  in  15  format. 
(Usually  this  will  be  body  1  but  it  can  be  any  body 

of  the  system. ) 

6.  Body  Connection  Array.  Next,  the  body  connection  array 
JFV(I)  (See  the  previous  section  for  a  description.)  is 
read  with  10  numbers  per  card  in  1015  format. 

7.  Acceleration  of  Gravity.  The  next  card  contains  the 
acceleration  of  gravity  to  be  used  by  the  program,  read 
in  F10.9  format.  The  choice  of  this  value  sets  the 
units  of  length  and  time.  All  other  input  data  must 
correspond  to  these  units.  (Since  the  dynamics 
formulation  is  independent  of  units,  the  output  will 
reflect  these  units  as  well.) 

8.  Body  Types  for  Masses  and  Inertia  Matrices.  The  next 
series  of  cards  is  used  to  read  or  input  the  masses  and 
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inertia  matrices  for  the  bodies  of  the  system.  In  many 
cases,  a  number  of  the  bodies  of  the  chain  system  are 
identical.  Hence,  the  first  card  in  this  series  contains 
a  number  KT  read  in  15  format  specifying  the  number  of 
different  types  of  bodies.  (If  KT=0,  the  bodies  are 
considered  to  be  all  distinct  for  the  purposes  of  input 
data,  even  though  some  may  be  physically  the  same.) 

9.  Masses.  If  KT=1,  all  bodies  are  considered  identical, 
each  having  the  same  mass.  In  this  case,  the  next  card 
contains  one  number  read  in  F10.9  format  signifying 
that  mass. 

If  KTj*l,  all  bodies  are  considered  to  have  different 
masses  (even  though  some  may  be  the  same) .  Hence,  in 
this  case,  the  masses  for  the  NB  bodies  are  read 
sequentially  in  the  next  series  of  cards  in  5F10.9  format. 

2 

Together  with  the  acceleration  of  gravity  (L/T  )  the 
masses  (M)  complete  the  choice  of  units  for  the  program. 
All  input  data  to  follow  must  be  consistent  with  this 
choice.  Note  also  that  the  output  of  the  program  will 
reflect  the  units  of  the  input  data. 

0.  Inertia  Matrices.  If  KT/0,  that  is,  if  there  are  several 
(at  least  1)  types  of  identical  bodies,  then  the  next 
card  contains  the  number  of  bodies  of  the  first  type 


followed  (on  the  same  card)  by  the  body  numbers  for  the 
bodies  of  the  first  type,  read  in  1015  format.  This 
card  is  immediately  followed  by  three  cards,  three 
values  per  card  in  3F10.9  format  specifying  the  inertia 
matrix  for  the  first  type.  This  procedure  is  then 
repeated  for  the  next  type  and  continued  until  all  KT 
types  are  exhausted. 

If  KT=0,  the  inertia  matrices  for  the  NB  bodies  of  the 
system  are  considered  to  be  distinct  (even  though  some 
may  be  the  same) .  In  this  case,  the  inertia  matrices 
for  the  NB  bodies  are  read  in  sequence  with  three  values 
per  card  in  3F10.9  format. 

11.  Mass  Center  Positions:  r  Vectors.  The  next  series  of 
cards  is  used  to  read  or  input  the  mass  center  positions 
for  the  bodies  of  the  system  —  the  r  vectors  (See 
previous  section) .  The  first  card  in  this  series  contains 
a  number  K  specifying  the  number  of  different  kinds  or 
types  of  identical  bodies,  with  respect  to  mass  center 
location,  read  in  15  format. 

If  K>0,  that  is,  if  there  are  several  (at  least  1)  types 
of  identical  bodies,  then  the  next  card  contains  the 
number  of  bodies  of  the  first  type  followed  (on  the  same 
card)  by  the  body  numbers  for  the  bodies  of  the  first 
type,  read  in  1015  format.  This  card  is  immediately 
followed  by  a  card  specifying  the  values  of  the  three 
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components  of  the  r  vector,  three  values  per  card  in 
3F10.9  format.  This  procedure  is  then  repeated  for  the 
next  type  and  continued  until  all  K  types  are  exhausted. 

If  K=0 ,  the  mass  center  positions  for  the  NB  bodies  of 
the  system  are  all  considered  to  be  distinct  (even 
though  some  may  be  the  same) .  In  this  case,  the  r 
vectors  are  read  in  sequence  with  three  values  (the 
components)  per  card  in  3F10.9  format. 

12.  Joint  Positions:  £  Vectors.  Similarly,  the  next 

series  of  cards  is  used  to  read  or  input  the  joint  (or 
reference  point)  positions  for  the  bodies  of  the  system 
(except,  of  course,  B^) .  These  are  the  NB-1  £  vectors. 
(See  previous  section.)  The  first  card  in  this  series 
contains  a  number  K  specifying  the  number  of  different 
kinds  or  types  of  identical  bodies,  that  is,  bodies 
with  identical  £  vectors,  read  in  15  format. 

If  K>0,  that  is,  if  there  are  several  (at  least  1) 
types  of  identical  bodies,  then  the  next  card  contains 
the  number  of  bodies  of  the  first  type  followed  (on  the 
same  card)  by  the  body  numbers  for  the  bodies  of  the 
first  type,  read  in  1015  format.  This  card  is  immediately 
followed  by  a  card  specifying  the  values  of  the  three 
components  of  the  £  vector,  three  values  per  card  in 
3F10.9  format.  This  process  is  then  repeated  for  the 
next  type  and  continued  until  all  K  types  are  exhausted. 
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If  K=0,  the  joint  positions  are  all  considered  to  be 
distinct  (even  though  some  may  be  the  same) .  In  this 
case  the  NB-1  £  vectors  are  read  in  sequence  with  three 
values  (the  components)  per  card  in  3F10.9  format. 

Magnetic  Tape  Card.  A  zero  or  blank  card  should  be  read 
here  if  no  magnetic  tape  is  to  be  used  to  record  the  output 
data.  If  a  magnetic  tape  is  to  be  used,  a  value  K  is  read 
in  15  format  indicating  the  output  unit  number. 

Heading  Card.  On  this  card  anything  may  be  written  as 
a  heading  for  the  particular  data  used.  This  heading 
appears  at  the  beginning  of  each  output  data  page. 

Specified  Motion  for  the  Bodies  of  the  System.  (See 
previous  section.)  For  the  acceleration  profiles  described 
in  the  previous  section,  the  data  cards  are  coded  as 
follows:  The  first  card  contains  the  number  of  profiles 
(up  to  10)  that  are  used,  read  in  15  format. 

The  next  cards  contain  the  number  of  points  (up  to  25) 
and  the  variable  numbers  of  the  respective  acceleration 
profile  curve  for  each  acceleration  profile  used.  They 
are  read  sequentially  in  1015  format. 

The  next  group  of  cards  contain  the  time  and  acceleration 
ordinate  for  each  point  of  the  respective  profile  curves 
which  are  used.  The  first  card  of  each  subgroup  (i.e.  for 
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each  profile  used)  also  contains  the  initial  velocity 
and  displacement  for  that  profile.  The  data  on  these 
cards  are  read  in  5F10.9  format. 

Blank  Cards.  Two  blank  cards  are  included  here. 

Gravity  Forces.  If  the  bodies  of  the  chain  system  are 
in  a  vertical  (x ^ ‘  gravity  field,  a  weight  force  is 
automatically  introduced  through  the  mass  center  of  each 
of  the  bodies  by  reading  1  on  the  next  card  in  15  format. 

If  gravity  iur  'es  are  not  to  be  included  automatically, 
a  blank  card  or  zero  is  read. 

Normal  Fluid  Drag.  If  the  cable  segments  of  the  chain 
system  are  immersed  in  a  fluid  bath,  a  normal  drag  force 
is  automatically  introduced  through  the  mass  center  of 
each  of  the  segments  and  a  normal  drag  torque 
is  automatically  applied  to  each  of  the  segments 
by  reading  a  1  on  the  next  card  in  15  format.  If  normal 
drag  forces  and  torques  are  not  to  be  included  automatically 
a  blank  card  or  zero  is  read. 

Tangential  Fluid  Drag.  If  the  cable  segments  of  the  chain 
system  are  immersed  in  a  fluid  bath,  a  tangential  drag 
force  is  automatically  introduced  through  the  mass  center 
of  each  of  the  segments  by  reading  a  1  on  the  next  card 
in  15  format.  If  tangential  drag  forces  are  not  to  be 
included  automatically,  a  blank  card  or  zero  is  read. 


20.  Added  Mass  Effect.  If  the  cable  segments  of  the  chain 
system  are  accelerating  in  a  fluid  bath,  an  added  mass 
force  is  automatically  introduced  through  the  mass  center 
of  each  of  the  segments  and  an  added  mass  torque 

is  automatically  applied  to  each  of  the  segments 
by  reading  a  1  on  the  next  card  in  15  format.  If  the 
added  mass  forces  and  torques  are  not  to  be  included 
automatically,  a  blank  card  or  a  zero  is  read. 

21.  Buoyancy  Force.  Buoyance  forces  are  automatically  intro¬ 
duced  through  the  mass  center  of  each  of  the  bodies  (cable 
segments  and  spherical  bodies)  by  reading  a  1  on  the  next 
card  in  15  format.  If  the  buoyancy  forces  are  not  to  be 
included  automatically,  a  blank  card  or  a  zero  is  read. 

22.  Two  Dimensional  Motion.  The  next  card  contains  a  dimension 
specification  switch,  read  in  15  format.  If  the  value 

1  is  read,  a  two-dimensional  dynamics  analysis  ensues. 

This  cuts  computer  run  time  required  for  the  problem 
solution,  but  is  otherwise  transparent  to  the  user. 

Motion  is  assumed  to  take  place  in  the  Z-Y  plane  and  the 
input  data  must  correspond  accordingly.  If  the  value  0 
(zero)  is  read,  the  standard  three-dimensional  analysis 
is  performed. 

23.  Known  and  Unknown  Variables.  In  most  problems,  some  of 
the  variables  will  be  specified  (known)  and  the  remaining 

1 9 


variables  are  to  be  determined  (unknown) .  To  identify 
these  variables,  proceed  as  follows:  Let  KT  be  the  number 
of  known  variables;  or  alternatively,  let  -KT  be  the 
number  of  unknown  variables.  Hence,  on  a  card  in  15 
format  read  either  KT  (positive  number  representing  the 
number  of  known  variables)  or  -KT  (negative  number 
representing  the  number  of  unknown  variables) . 

Next,  in  1015  format  using  as  many  cards  as  needed,  read 
the  variable  numbers  (e.g.  7  for  X?,  etc.)  for  the  KT 
known  variables  or  the  -KT  unknown  variables,  according 
to  whether  KT  or  -KT  is  read  in  the  previous  card. 

Finally,  if  KT=0  (i.e.  a  blank  card) ,  all  variables  are 
considered  unknown. 

24.  Initial  Values  of  the  Variables.  The  initial  values  of 

the  variables  are  included  as  input  whether  the  variables 
are  known  (specified)  or  unknown  (to  be  determined) . 

If  the  variables  are  known  and  the  initial  values  are 
different  than  the  initial  values  of  the  specified 
function,  say  in  SUBROUTING  START,  then  the  latter  will 
be  given  preference  in  the  program. 

It  is  convenient  to  consider  the  variables  and  hence,  their 
initial  values  in  groups  of  three,  with  the  first  group 
being  the  displacement  of  in  reference  frame  R.  This 
group  is  designated  as  the  "R-group"  or  "body-zero"  group. 
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The  ordering  of  the  variables  within  the  group  is  the 
natural  or  usual  order:  x^=x,  X2=y,  and  x^z.  The  next 
group  of  three  variables  are  the  orientation  angles  of 
Body  1  and  they  are  designated  as  the  group.  The 
ordering  of  the  variables  within  the  group  is  again  the 
natural  or  usual  order  with  x4»a^,  x5=®i»  and  X6~Y1' 

This  procedure  is  continued  through  the  3NB+3  variables 
and  the  NB  bodies,  so  that  in  general  the  group  of 
variables  associated  with  body  B;  (the  "B  ^ "  group)  are  the 
three  variables  defining  the  orientation  of  B ^  relative 
to  its  adjacent  lower-numbered  body,  with  the  usual 
ordering:  a,S,Y- 

There  are  several  options  available  for  the  input  of  the 
initial  values  of  these  variables:  i)  They  may  be  read 
in  sequence;  ii)  They  may  be  read  at  random  (in  groups  of 
three) ;  or  iii)  They  may  be  read  in  groups  of  three  — 
when  identical  values  are  involved. 

To  actuate  one  of  these  options,,  read  a  number  K  in 
15  format  on  a  card,  with  the  input  data  following  on 
subsequent  cards  according  to  the  following  direction: 

If  K*0  (or  a  blank  card) ,  the  initial  values  are  read 
sequentially  in  groups  of  three  with  three  values  per 
card  in  3F10.9  format  (from  "body-zero"  to  "body-NB" 

(NB+1  cards)). 
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If  K<0,  the  initial  values  of  the  variables  are  read  in 
-K  groups  of  three,  but  the  groups  themselves  may  occur 
in  random  order.  Also,  if  fewer  than  the  NB+1  groups  are 
read,  the  initial  values  of  the  variables  associated  with 
the  groups  not  read  are  automatically  set  equal  to  zero. 

To  identify  the  groups,  the  first  value  on  each  card  is  the 
"body  number"  of  the  group  in  15  format.  The  initial  values 
of  the  variables  are  then  read  in  F15.9,  2F10.9  format  on 
the  same  card  with  a  total  of  four  values  per  card. 

If  K>0,  the  initial  values  of  the  variables  are  read  in 
K  groups  .of  groups  of  three.  (This  option  is  used  when  the 
initial  values  of  the  variables  between  two  or  more  groups, 
or  bodies,  are  identical.)  The  cards  are  coded  as  follows: 
On  the  first  card  (following  the  card  with  K) ,  read  the 
number  of  groups  (or  bodies,  with  identical  initial  values 
of  the  variables)  within  the  first  group.  Then  read,  on 
the  same  card,  the  corresponding  body  numbers.  All  numbers 
on  this  card  are  in  1015  format.  This  card  is  immediately 
followed  by  a  card  specifying  the  three  initial  values 
of  the  variables  in  3F10.9  format.  This  procedure  is 
then  repeated  for  the  next  group  and  continued  until  all 
K  groups  are  exhausted. 

The  second  option  above  is  particularly  helpful  for  making 
the  initial  values  of  all  the  variables  zero.  To  do  this, 
simply  read  K  as  -1  and  follow  this  with  a  card  having  a 
value  1,  read  in  15  format. 


25.  Initial  Value  of  the  Derivatives  of  the  Variables.  These 
are  read  in  exactly  the  same  manner  as  the  initial  values 
of  the  variables  as  described  in  24.  above. 

26.  Printing  Priority.  On  a  single  card  in  15  format  a  number 
is  read  between  -100  and  +100  determining  the  amount  of 
print-out  for  a  given  run  as  follows : 

If  0,  all  print-out  (the  "standard  print-out")  for  every 
printing  time  (see  27.  below)  is  given. 

If  >10  printing  of  X  and  DX  in  radians  is  omitted. 

If  >20  printing  of  moments  between  bodies,  FF  and  all 
fluid  induced  forces  are  omitted. 

If  >30  printing  of  DDX  in  degrees  is  omitted. 

If  >40  printing  absolute  and  relative  velocity  and 
acceleration  is  omitted. 

If  >50  printing  of  mass  center  and  joint  position  is  omitted. 
If  >90  printing  of  X  and  DX  in  degrees  is  omitted. 

If  >95  All  printing  is  omitted. 

If  >0  (negative)  the  corresponding  print-out  is  given  for 
all  integration  steps,  and  all  print-out  at  every  print 
time  (see  27.  below)  is  given.  (For  example,  -100  gives 
the  maximum  possible  print-out.) 
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27.  Printing  Time.  This  is  a  value  read  on  a  single  card  in 
F10.9  format  determining  the  time  increment  for  print-out. 
This  may  be  the  same  or  different  (usually  greater)  than 
the  integration  time  increment  PRMT(3). 

28.  PRMT.  This  is  an  array  of  4  numbers  needed  for  the 
integration  subrouting  RKGS.  They  are  entered  on  a  single 
card  in  5F10.9  format.  The  first  number  is  the  integration 
starting  time,  the  second  number  is  the  integration  ending 
time,  the  third  number  is  the  integration  time  interval, 
and  the  fourth  number  is  the  desired  upper  bound  on  the 
error  of  the  integration. 

If  the  third  number  PRMT (3)  is  read  as  zero,  then  PRMT (3) 

is  automatically  set  to  0.1.  If  PRMT (3)  is  read  as  negative, 

— M 

say  -M,  then  PRMT (3)  is  automatically  set  to  2 

If  PRMT (4)  is  read  as  zero,  then  PRMT (4)  is  automatically 
set  to  0.001. 

29.  Fluid  Data.  The  next  three  cards  are  used  to  read  fluid 
related  data.  The  first  card  contains  the  inertial 
components  of  the  fluid  velocity  in  a  5F10.9  format.  The 
second  card  contains  the  inertial  components  of  the  fluid 
acceleration  in  a  5F10.9  format.  (Note:  Since  the  fluid 
velocity  is  currently  assumed  constant  in  space  and  time,  the 
components  of  the  fluid  acceleration  must  be  zero.)  The 
third  card  contains  values  for  the  Y  coordinate  of  the 
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surface  height  of  the  fluid,  the  density  of  the  fluid, 
and  the  viscosity  of  the  fluid  in  5F10.9  format. 


30.  Blank  Card.  A  blank  card  is  included  here. 

31.  Number  of  Towed  Spheres  and  Spherical  Anchors.  The  next 
card  contains  the  number  of  spherical  bodies  NSPHRB,  read 
in  15  format.  Note:  These  bodies  may  only  be  at  the  end 
of  a  branch. 

32.  Body  Numbers  of  Towed  Spheres  and  Spherical  Anchors. 

The  following  card(s)  contain  the  body  number (s)  of  the 
spherical  bodies,  read  in  1015  format.  If  NSPHRB  the 
number  of  spherical  bodies  is  2ero  these  cards  are 
omitted . 

33.  Body  Data.  The  next  Nb  cards  contain  body  data  for  each 
of  the  bodies  of  the  chain  system.  There  is  one  card  for 
each  body  in  sequence  from  body  1  to  body  NB.  If  the 
body  is  a  cable  segment  then  its  data  card  contains  its 
length,  its  diameter,  and  the  segment  fixed  components  of 
a  unit  vector  parallel  to  the  centerline  of  the  segment, 
read  in  5F10.9  format.  (Note:  The  cable  segments  are 
assumed  to  be  cylindrical  in  shape.)  If  the  body  is  a 
spherical  body,  then  its  data  card  contains  its  diameter, 
read  in  5F10.9  format. 
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34.  END  OF  DATA  Card.  The  last  card  must  read  "*END  OF  DATA*". 


(See  example  input  data.) 


V.  EXAMPLE  LISTING  OF  INPUT  DATA 

Figure  4.  contains  an  illustration  of  a  cable  model.  An 
example  listing  of  the  input  data  for  this  particular  model  is 
contained  in  the  following  pages. 
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CONTENTS  OF  CARD 


VI .  OUTPUT  DATA 


The  output  data  is  labeled  on  the  computer  print-out  and 
thus  is  primarily  self-explanatory.  It  consists  of  two  parts: 

The  first  part  is  simply  an  "echo"  or  copy  of  the  input  data. 

(In  the  input  motion  section,  the  velocity  and  position  at  the 
beginning  of  each  acceleration  interval  is  computed  and  also 
listed.)  The  second  part  contains  at  each  output  time  (depending 
upon  the  printing  priority  (See  26,  in  the  input  section.)): 

i)  All  variables  and  the  first  derivative  of  the  variables 
(angles  are  in  radians) . 

ii)  Joint  positions  (in  both  inertia  space  and  relative  to 
the  reference  body) . 

iii)  Mass  center  positions  (in  both  inertia  space  and  relative 
to  the  reference  body) . 

iv)  Mass  center  velocities  and  accelerations. 

v)  All  variables  and  the  first  and  second  derivatives  of 
the  variables  (angles  are  in  degrees) . 

vi)  Moments  and  forces  associated  with  variables  which  are 
specified. 

vii)  Fluid  induced  forces  and  torques . 

An  example  listing  of  output  is  contained  in  the  following 
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VII.  COMMENTS  ON  UNITS  AND  ORIENTATION  ANGLES 


The  units  most  frequently  used  in  the  program  are: 

1)  slug-in-sec;  2)  slug-ft-sec;  and  3)  kg-m-sec.  However,  any 
other  set  of  consistent  units  may  be  used.  The  set  of  units 
the  program  uses  is  determined  strictly  by  the  input  data. 
Whatever  set  is  chosen  by  the  user  must  be  used  for  the  input 
data,  then  the  program  will  use  this  set  as  well. 

There  is  a  singularity  present  in  orientation  angles.  The 
6  rotation  angle  should  not  be  allowed  to  go  through  90°.  This 
can  usually  be  accomplished  by  properly  defining  the  initial 
position  of  the  model  so  that  explicit  motions  relative  to  this 
initial  position  are  such  that  6  does  not  go  through  90°. 


36 


4 


APPENDIX 


ORIENTATION  ANGLES 


As  mentioned  in  the  foregoing  text,  a  majority  of  the  depen¬ 
dent  variables  are  angles  measuring  the  relative  orientation  of 
adjacent  bodies  of  the  chain  system.  These  angles  are  defined 
as  follows:  Consider  two  coincident  dextral  axes  systems  X^,  Y^, 

Z^  and  X2 ,  Y3/  Z2  as  shown  in  Figure  6.  Let  X^,  ,  Z^,  be  fixed 

in  one  of  the  adjacent  bodies  (the  lower  numbered  body) ,  and  let 
X2>  Y2/  Z2  be  free  to  rotate  relative  to  X^,  Y^  with  the 
X^  and  X2  axes  remaining  coincident.  The  angle  a  is  then  defined 
to  be  the  dextral  rotation  of  Y 2 ,  Z2  relative  to  X^,  Y^,  Z^ 

as  shown  in  Figure  7.  Next,  let  X^,  Y3,  be  a  third  dextral 
axes  system  originally  coincident  with  X2,  Y2,  Z2*  Let  X3,  Y3,  Z3 
rotate  relative  to  X2,  Y2»  Z2  with  Y2  and  Y3  remaining  coincident. 
The  dextral  rotation  of  X3,  Y3,  Z3  about  Y3  and  Y3  defines  the 
angle  8  as  shown  in  Figure  8.  Finally,  let  X4 ,  Y^ ,  Z4  be  a 
fourth  dextral  axis  system  originally  coincident  with  X3,  Y3,  Z3 
and  free  to  rotate  relative  to  X3,  Y3,  Z3  with  the  Z3  and  Z4  axes 
remaining  coincident.  A  dextral  rotation  about  Z3  and  Z4  defines 
the  angle  y  shown  in  Figure  9.  Let  X4 ,  Y4 ,  Z4  be  fixed  in  the 
other  adjacent  body  (the  higher  numbered  body) .  The  three  angles 
a,  8,  and  y  thus  measure  the  relative  orientation  of  the  two  bodies. 

In  the  foregoing  text  the  variables  are  listed  in  groups  of 
three.  Those  variables  which  are  orientation  angles  are  sequen¬ 
tially  "a",  "8"  and  "y"  angles  as  defined  above. 
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20.  AMTRACT  (CmtHmm  m  ro*m—  o!4m  It  iiMiiiay  m*4  ibomtlty  by  M  mA  nwMr;  Vhi.S  report  presents  a 

computer  program,  UCIN-CABLE:  Underwater  Cables.  The  program  is  de 
signed  to  study  the  dynamics  of  underwater  cables.  The  cable  is 
modeled  as  a  series  of  rigid  cylinders  connected  end-to-end  by  ball 
and  socket  joints.  There  may  be  one  or  many  branches,  and  no  close 
loops  are  formed.  The  effects  of  normal  fluid  drag,  tangential  -lu 
drag,  added  mass,  and  buoyancy  are  included.  This  manual  provides 
detailed  instructions  for  building  a  set  of  input  data  and  inter- 
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